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Abstract. 

We consider N particles interacting pair-wise by an inverse square potential in 
one dimension (Calogero-Sutherland-Moser model). When trapped harmonically, its 
classical canonical partition function for the repulsive regime is known in the literature. 
We start by presenting a concise re-derivation of this result. The equation of state is 
then calculated both for the trapped and the homogeneous gas. Finally, the classical 
limit of Wu's distribution function for fractional exclusion statistics is obtained and we 
re-derive the classical virial expansion of the homogeneous gas using this distribution 
function. 
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1. Introduction 

We consider a system of identical particles in one dimension interacting via an inverse 
square pairwise interaction. This is a class of integrable many-body systems known as 
the Calogero-Sutherland-Moser (CSM) model [H |2] which is a classic example of an 
exactly solvable many-body system. For over three decades, the variants of this model 
have provided a template for analysing disparate problems in high energy and condensed 
matter physics. 

The Hamiltonian of the model in the presence of a harmonic confinement is given 

by 

where Xi^pi denote the positions and momenta of the particles, and A is dimensionless 
coupling constant. 

In the absence of a harmonic confinement, the system is classically integrable [2]. 
Defining Ti^X = a^, the integrals of motion are constructed using the Lax matrix defined 

as 

Lij = PiSij + (1 - Sij) _ . (2) 

The integrals of motion are given by 

In = tr{L''), for n = l,...,N. (3) 

It is straightforward to show that the integrals of motion are in involution, [/„, /„] = 0, 
and hence the system is classically integrable. 

In this paper we are interested in the classical statistical properties of a system of 
identical particles whose dynamical behaviour is described by the Hamiltonian in Eq. 
([T]). The quantum dynamics of such particles has been extensively studied [HE], and its 
exact N-particle quantum canonical partition function is known. By taking its ^ 
limit, the corresponding classical partition function was derived long back [1]. More 
recently, further studies associated with the classical integrals and related Jacobians 
have been made O [6]. In the present paper, the emphasis is on the thermodynamic 
properties of this classical system, particularly in relation to the fractional exclusion 
statistics (FES) that it obeys in the quantum regime [3], [7]. In sect. 2, we give a concise 
derivation of the classical N-particle canonical partition function, and the corresponding 
equation of state for the trapped gas. The harmonically confined system has a constant 
density of states, and only the second virial coefficient is found to be non-zero. But 
our primary interest is in obtaining the virial expansion of the unconfined gas in the 
thermodynamic limit. The classical partition function for such a system is divergent. 
In sect. 3, we use the harmonic oscillator confinement as a regulator, and obtain the 
virial expansion for the equation of state in the limit of the oscillator frequency going to 
zero. In sect. 4, in the context of FES, we obtain the non-trivial classical distribution 
function by taking the appropriate limit of Wu's quantum occupancy factor for particles 
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obeying fractional exclusion statistics [11 E]- Using this, we derive the virial expansion 
of the homogeneous gas and confirm the results obtained in sect. 4. We also find the 
energy of the classical CSM model at zero temperature. 



2. Classical limit of the quantum canonical partition function 

The classical partition function for identical particles is given by 

Zn^ = j^^J^^^N J d""^ d'^p expi-PHr,), (4) 

where the Hamiltonian is given by Eq. ([1]). We choose the interaction strength 
h^X = o? ^ with a real. The interaction is therefore always repulsive. Note the explicit 
factor of iV! in the above expression for identical particles. While it is rather difficult to 
obtain the classical partition function by direct integration for all A^, we may approach 
the problem as the classical limit of the quantum canonical partition function Zm{I3). 
This was the route that was taken originally in [Ij. We obtain the desired result more 
directly, however, by using a property of CSM that relates to FES. To this end, we set 
the interaction strength A = g{g — 1), where > 0, and recall the known result [9] 

Z, = e^^-(-^)^Z-, (5) 

where Zj^ is for an arbitrary g, and Zfj is the A^-particle non-interacting fermion 
partition function (for = 1). Identical results are obtained when the problem is 
approached from the bosonic representation (for g = Q). 

We are now in a position to take the classical limit. We define the classical limit 
as one in which the parameters m,uj,f],a are held fixed and we take the limit h 0. 
Assuming this we proceed as follows. Note that we have set 

= h'X = h'g{g - 1), (6) 

where the classical interaction strength a is independent of h. Both a and g are positive 
definite and are related by 

h 



gn=- + ^a' + h'/A. (7) 
Now taking the limit h {] while keeping a fixed implies g ^ oo and 

gTi a. (8) 

This defines the approach to the classical limit. Taking this limit in Eq. ([5]) we obtain 

limZ^v = e-"^"^^^ lim[Zi^]. (9) 



N(N-\) 
' 2 

The non-interacting fermion partition function is given by the well-known 
expansion, namely. 



V I 



I 



^ - ^^^^z,mzr\P) + ■■■], (10) 
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where the sum over V is given by the number of partitions of such that E/li "^i^ = ^ 
for ui and / positive integers and 

^'(« = 2sinh(L/2) - 

Now, taking the hmit — 0, we obtain the desired result 

a JV(iV-l) 1 

This is the same resuh obtained in ^ using a different but longer method. 
2.1. Equation of state of the trapped classical gas 

We can easily calculate the thermodynamic properties of the trapped gas from Zm{(3) 
given in Eq. (IT^ . The free energy is given by F/v = — rlnZ^v, where r = 1/13. Since 
the density of states is constant for harmonic confinement, it is like a two-dimensional 
gas, and the pressure is given by 

where A = = h/muj, I being the oscillator length. Writing 

(OfA ^ (mA (duA 

[dA)^ [du)^ [dA)' 
and p = N/A, we obtain the equation of state 

l3P = p + -^p\ (15) 
where P is the pressure. Note that virial coefficients of order three and higher are zero. 



(13) 



(14) 



3. Classical Equation of state for the homogeneous gas {u 0) 

We now consider the virial expansion of the unconfined gas in the thermodynamic limit. 
The classical partition function for such a system is divergent. However, we use the 
harmonic oscillator confinement as a regulator to calculate the partition function as 
before and obtain the virial expansion for the equation of state in the limit of the 
oscillator frequency going to zero. 

In the dilute limit, the equation of state is given by 

oo 

PP = Y1 ^'-^^ wliere z = e^'', (16) 
1=1 

and bi are the cluster expansion coefficients [lOl |TT] which appear in the fugacity 
expansion as above. They can be expressed in terms of the canonical partition functions 

bi = {z,y-' E (-i)(^-'"'-^nE^^ - 1)' n 

{nii} i i 



z\ 



(17) 
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Note that a cluster coefficient of order / involves all the partition functions up to Zi. The 
summation over mj is constrained by the number of partitions of /, that is Y^\=i infii = I. 
Using the expansion for the density in terms of the cluster coefficients, namely, 

oo 

p = Y.lhizK (18) 
1=1 

along with Eq. ( JT6l) . the virial expansion coefficients are defined by 



/3P = p [1 + ^a„(ATp)"-^], (19) 



where = '^2Tih^l3/m is the thermal wavelength. The are the virial coefficients of 
the system. Note that 6i = Oi = 1. 

Consider the limit in which the confinement is removed. We do this by taking the 
limit cij — > as follows. The cluster coefficients are given by. 

Note that the numerical pre-factors in these expressions for the harmonic regularisation 
— > are different from the box regularisation L — oo. In c? dimensions (rf = 1 in 
our case), the expressions for 6„ in the simple harmonic regularisation must be taken 
to be larger by a factor of ra'^/^, where n is the order of the virial coefficient, than for 
box regularisation. This ensures that they give the same result as — and L — oo 
respectively [12j . 

Then the first two virial coefficients are given by 

as = Ahl - 263. (21) 
For the quantum gas, we find that 

and therefore 

«2 = -^{g - 1/2), 

as = (2 - V^)g{g-l) + 1 - -1=. (23) 

For g = (1), we recover the virial coefficients for a one- dimensional gas of non- 
interacting bosons (fermions). If we take the limit ^ — > and g —>■ 00 keeping a 
fixed as before, we get the virial expansion for the classical gas 



PP = P[l + ^ ap + (2 - v^) ^ aV + ...]. (24) 



m m 
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4. Classical distribution function 

An interacting system of particles described by CSM model may be mapped on to an 
ideal gas obeying fractional exclusion statistics (FES) [3 [3]. In this section, we derive the 
classical virial expansion of the homogeneous gas from this starting point and confirm 
that we obtain the same results as in sect. 4. It is also shown that the energy per unit 
length at zero temperature may be consistently obtained from our classical description 
by taking the limit % 

The definition of the statistical parameter in FES, denoted by g{> 0), is based on 
the rate at which the number of available states in a system of fixed size decreases as more 
and more particles are added to it. The statistical parameter g assumes the values and 
1 for bosons and fermions respectively, because the addition of one particle reduces the 
number of available states by g. The application of the finite temperature distribution 
function [8j then enables us to calculate the temperature dependent quantities of the 
system. 

As is well known, the Haldane-Wu statistics is realized by the CSM model in 
one dimension [3], with the statistical parameter g in the FES being identical to the 
interaction strength in the CSM model as noted earlier. The potential and kinetic 
energy scale in the same way in this model, and both the energy densities scale as p^. 
The distribution function or average occupancy for FES particles has been derived by 
Wu |8] and is given by 

n, = (25) 
Wp + g 

where p denotes the momentum, and the dispersion relation is given by = /2m. 
The parameter g is called the statistical parameter of FES since the occupancy of a 
given momentum state depends on g. It has been shown that the statistical parameter 
of FES is also the interaction coupling in the CSM model as used in the previous section. 
The function Wp satisfies the equation 

wl{l + Wpf-' = e^^'^-^'\ (26) 

We now ask the question, what is the classical limit of the distribution function or, 
equivalently, how do we take the limit g oo? Let us assume that 



w. 



p 



ghp and e'^^ = g e-~^^% (27) 



where 7p depends on the momentum p, and /ic is a renormalised chemical potential 
relevant to the classical limit. They also depend on other variables like temperature 
and density as will become clear below. 

Using the above, we can write Eq. (l26l) in the following form 
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If we now take the limit g ^ oo keeping all other variables /3, 7^, e^, /ic fixed, we find 
that 

- — = e^(^^-^=\ (29) 
7p 

where we have made use of the identity linig^^il — ^pjgY = e~'^^. Note that Eq. 
( !29|) uniquely fixes 7^ in terms of the momentum p, temperature and the chemical 
potential since the function monotonically goes from 00 to as a function of 7^ for 
< 7p < 00. Furthermore, using Eq. (12^ and Eq. ^17) we see that 

which is now the analogous classical distribution function for CSM particles. 

There is yet another way of obtaining the above result from the grand partition 
function corresponding to FES, namely, 

lnZG = El^(l + -) = E- (31) 
p p 9 

in the limit g ^ 00. We immediately find that 

I dlnZc _l 7p 

[3 dep fi- 1 + 7p ' 

as desired. 

We can now obtain the classical limit by setting gh = a. In a system with length 
L, the free energy F = —{1/(3) InZc follows from Eq. (!3T|) : 

L f°° dp 'jp L f°° dp 



= —o^^ = ZT^^ (32) 



(3 J-00 2Tih g 13 J -00 2'Ka ^^^^ 
Since the pressure is given by P = —{dF/dL)j3, we see that 

The density (i.e., the number of particles per unit length) is given by 

dp f°° dp 7„ 

— T n-r, = 
2Txh ^ I 

while the energy per unit length is 

E=r^^]L. (36) 
27ra; 1 + 7p 2m 



00 



2^ "- = 7-00 2;^ iT^' ^''^ 



The virial expansion at high temperature can be obtained using Eqs. fl29|) . fl34j) and 
fl5S]) as follows. We find that as /5 — > 0, we must take e^^" ^ so that e-^^'-v-i^<^) < 1 
for all values of p. Using Eq. (129|) . we can expand 7p as a power series in e~^^''p~^''\ To 
go up to the third virial coefficient, we find that 

Q 
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Eq. (1531) then gives 



p = i I _ ^e^P^^^ + ^e=^^^^ + ■ • • V (3^ 



This equation can be inverted to give 



e 



- , , 72 J % (4 _ 5^) f . I , . . . .(39) 



Eqs. ( l34l) and ( 1371) now give 



^ f e/^A^= - ^e2/5^= + ^e^'^'^^ + ■ ■ ■ V (40) 
u: y 27113 \ v2 2 y 

Substituting Eq. ( l39l) in ( l40l) . we obtain the expression in Eq. (|2^ . 

Finally, let us consider the zero temperature limit. Note that as /3 — > cxd, 7^ = if 
Cp > pc and 7p = cx) if ep < /Xc; thus 7p/(l + 7p) is or 1 in these two cases. This is 
very similar to the Fermi distribution function at zero temperature. Using this fact in 
Eq. (|25l) . we find that there exists a Fermi momentum pp = yJ2mpc which is related to 
the density through p = pp/{na). Eq. (15^ then shows that the energy per unit length 
is given by = 7r^a;^p'^/(6m). Let us now show directly that this is the expected value 
of the classical energy at zero temperature. At T = 0, the particles are at rest; hence 
the kinetic energy is zero. The repulsive two-body interactions in Eq. ([1]) (where we 
have taken ^^A = as usual and also set lu = 0) will be minimised if the particles are 
equally spaced on a line, with the nearest neighbour spacing being equal to 1/p. If the 
particles are ordered such that Xj < Zj+i for all i, we will have — Xi = n/p. The 
interaction energy per particle is then given by 

_ y ! = 71 a p ^ ^^^^ 

m ~^ {n/pY 6m 

Thus the energy per unit length is given by 7r^Q;^p'^/(6m). 



5. Summary 

In this paper, we have used the exact solvability of the energy spectrum of the quantum 
CSM model for any value of the interaction parameter g to study the classical limit; 
this limit is obtained by taking ^ — >■ and g ^ 00 keeping gh = a fixed. Our 
derivation of Zi^{(3) is more concise than previous derivations. We have computed 
the virial expansions for the classical CSM model with or without a harmonic confining 
potential (i.e., for a homogeneous system). Finally, we have found the classical limit of 
the Wu distribution function for FES and used this to show consistency between the 
virial expansions of the homogeneous CSM model obtained from Z[^{(3) and from FES 
up to the third virial coefficient. 
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